Spiral graphone and one sided fluorographene nano-ribbons 
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The instabihty of a free-standing one sided hydrogenated/fluorinated graphene nano-ribbon, i.e. 
graphone/fluorographene, is studied using ab-initio, semiempirical and large scale molecular dy- 
namics simulations. Free standing semi-infinite arm-chair like hydrogenated/fluorinated graphene 
(AC-GO/ AC-GF) and boat like hydrogenated/fluorinated graphene (B-GO/B-GF) (nano-ribbons 
which are periodic along the zig-zag direction) are unstable and spontaneously transform into spiral 
structures. We find that rolled, spiral B-GO and B-GF are energetically more favorable than spiral 
AC-GO and AC-GF which is opposite to the double sided flat hydrogenated/fluorinated graphene, 
i.e. graphane/fluorographene. We found that the packed, spiral structures exhibit unexpected local- 
ized HOMO-LUMO at the edges with increasing energy gap during rolling. These rolled hydrocarbon 
structures are stable beyond room temperature up to at least T=1000K. 



I. INTRODUCTION 

The discovery of graphene has been a driving force 
for the scientific community to synthesize and to char- 
acterize new materials with similar morphologies due to 
their unique properties [5-^. Fully double sided hydro- 
genated graphene (GE), i.e. graphane (GA), and double 
sided fluorographene are quasi two-dimensional lattices 
of carbon (C) atoms ordered into a buckled honey- 
comb sublattice, where each carbon atom is covalently 
bonded to hydrogen (H) or fluor (F), respectively, in an 
alternating, chair-like arrangement [6, 7]. The chemisorp- 
tion of hydrogen and fluor atoms results in an important 
reconstruction of the chemical bonds and angles of the 
underlying honeycomb lattice [8]. This transition from 
sp^ to sp^ hybridization turns the graphitic C-C bonds 
into single bonds by the formation of additional single C- 
H and C-F bonds, which change locally the planar shape 
of graphene into an out of plane, angstrom scale buckled 
geometry [9]. 

Experimentally, it has been shown that GA can be ob- 
tained reversibly starting from a pure GE layer in the 
presence of atomic hydrogen [10]. It has become an in- 
teresting material due to its potential applications in na- 
noelectronics [11]. On the other hand, the experimentally 
obtained single-layer fluorographene exhibits a strong in- 
sulating behavior with a room temperature resistance 
larger than 1 Tl^, and a high temperature stability up 
to 400 [7]. 

Both in experiments [TO] and in ab-initio calculations, 
fiat single-sided hydrogenated graphene, graphone (GO), 
was found to be unstable [l2l - [l7| which was also 
demonstrated by phonon band structure calculations. 
Supporting graphene on a substrate and hydrogenating 
was demonstrated to produce GO in a recent experi- 
ment [l^. Ab-initio calculations showed that arm-chair 
(AC)-GA is more stable than boat like (B)-GA while 
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FIG. 1: (Color online) Side views of a small portion of simu- 
lated graphones which were relaxed at T=10K using molec- 
ular dynamics simulation. Arm-chair graphone (a) and boat- 
like graphone (b). 



B-GO was found to be more stable than AC-GO [19|. 
Ab-initio calculations also showed that GO exhibits mag- 
netism due to the localized electrons on the carbon atoms 
without hydrogens, in contrast to nonmagnetic graphene 
and graphane [20l, l2l| . 

By using both ab-initio calculation and elasticity the- 
ory, Kudin et al found that while linear elasticity theory 
can be used for studying the stiffness and flexural rigidity 
moduli of carbon and boron-nitride nanoshells, one-sided 
fluorinated carbon ribbons exhibit a strong tendency to 
shell formation resulting into very small diameter tubes 
with (4,4) and (3,3) indexes and properties that devi- 
ate from linear elasticity theory [22]. Such a strong ten- 
dency to bending is also observed for Si and bilayer SiGe 
nanofilms, due to reconstructions in the surface layer [2_3,] 
associated with dimer formation. Although the possible 
deformation of the edge of graphene (see [24] and refer- 
ences therein) and the synthesis of carbon nanotubes via 
one-sided hydrogenation or fluorination of supported car- 
bon films have been investigated theoretically [25^, the 
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geometrical structure of long (referring to the distance 
between the free edges) nano-ribbons of graphone and 
their corresponding stability and electronic properties at 
room temperature have not been studied so far. Here we 
show that such long nano-ribbons of GO/GF do not form 
carbon nanotubes, as in the case of relatively short sam- 
ples with specific sizes [25], but, instead, spontaneously 
form rolled, spiral structures with interesting localization 
of frontier molecular orbitals. The one-sided absorption 
and corresponding sp^ hybridization is the main driving 
force for this transformation, and the final spiral struc- 
ture is significantly more stable than the corresponding 
fiat structures. We used the three mentioned methods 
to investigate various electronic and geometrical proper- 
ties of these new spiral structures, which we named spiral 
graphone (fluorinated graphene)\ 

There have been several recent experimental and the- 
oretical studies related to the formation of carbon nan- 
otubes using a graphene sheet and inverse processing [26- 
[28j . Many experimental evidences point towards the 
possibility of the production of the here proposed spi- 
ral carbon nanostructures. In a recent review paper, Vi- 
latela and Eder |29| reviewed related experimental stud- 
ies, compared the potentials and characteristics of nano 
carbon composites (such as various rolled sheets) and 
discussed key challenges for the use of the new carbon 
nanostructures. Furthermore, similar rolled structures 
for Ni3Si205(OH)4 were synthesized and analyzed using 
X-ray diffraction (XRD) and transmission electron mi- 
croscopy (TEM) [30]. Shen et al [31] reported recently 
the synthesis of self- assembled kinked In20s nanospi- 
rals and multikinked nanowires using a tube-in-tube laser 
ablation chemical vapor deposition (CVD) method with 
gold nanoparticles as the catalysts. 

Apart from being more stable than the fiat structures, 
we also show that spiral B-GO and B-GF are more sta- 
ble than spiral AC-GO and AC-GF, respectively. The 
spiral structures are closely packed and stable even at 
1000 K and are candidates for a new class of carbon 
nano-structures. Since the deformations for the GFs 
and the GOs are very similar we will report here mainly 
the results for the GOs. After rolling, the highest oc- 
cupied orbital (HOMO) and the lowest unoccupied or- 
bital (LUMO) are separated and appear to be localized 
at opposite (free) ends, making these two ends chemically 
more active [32^ , while before and during rolling they are 
localized at both ends simultaneously. We found that the 
rolling process makes the systems strongly polarized and 
more insulating. 



II. METHODS AND MODELS 

We performed ab-initio calculations using GAUSSIAN 
(G09) [33|] which is an electronic structure package that 
uses a basis set of Gaussian type of orbitals. In the ab- 
inito calculations for the exchange and correlation (XC) 
functional, the hybrid functional B3LYP is adopted in 



G09. The self consistency loop was iterated untill the 
change in the total energy is less than 10~^ eV, and the 
geometries were considered relaxed once the force on each 
atom is less than 50 meV/A. Using the 6-31G* basis set in 
G09, we expect that our calculation is capable to provide 
a reliable description of the electronic properties of the 
different systems. In the case of AC-GO the A-sublattice 
sites of GE were covered with hydrogens, while the B-GO 
model was obtained by hydrogenating all carbon atoms 
involved in horizontal C-C bonds along the AC direction. 
Hence the number of carbon atoms is twice as large as 
the number of hydrogen atoms. The simulated samples 
(using the DFT method) have typically more than 500 
atoms. 

For the larger samples, we performed semi-empirical 
calculations (PM6 level of calculation in G09) as well as 
classical molecular dynamics simulations (MD) using the 
LAMMPS package [34]. 

For the classical MD simulations, the modified second 
generation of Brenner's bond-order potential, i.e. Adap- 
tive Inter molecular Reactive Bond Order (AIREBO) [35| 
and the ReaxFF [36] potentials were employed to simu- 
late the GOs and GFs, respectively. These simulations 
were performed using periodic boundary conditions along 
the lateral side (perpendicular to the rolling direction) in 
the canonical NVT ensemble with a Nose-Hoover ther- 
mostat for temperature control and a time step of 0.1 fs. 
The simulation box size along the rolling direction was 
taken equal to twice the system size (in its fiat geometry) 
in that direction. 

AIREBO consists of two parts, namely the reactive 
bond order potential (REBO [37]) for the short range in- 
teractions (<2 A), but with a modified, four body, bond 
order contribution for the torsion in various hydrocarbon 
configurations, and the van der Waals (vdW) term for 
long range interactions within distances 2 A < r < 10 A 
which is similar to the standard Lennard Jones potential, 
but with an environment dependent (adaptive) suppres- 
sion of the (too) strong 1/r-^^ repulsion. Its total energy 
reads: 

E AIREBO = Ere BO + ^vdW, (1) 

where the prime in E'j^^^q is added to indicate the mod- 
ification in the torsion term. In our MD simulations 
we checked the effect of the vdW term (by turning it 
on/off) to insure that when GO is rolled up the vdW term 
does not prevent possible bonding between two close edge 
atoms. The simulations where we included the vdW term 
are labeled with the symbol *. We found that including 
the vdW term lowers the energy of spiral GO* by only 
0.04eV/atom with respect to GO. 

In ReaxFF (for simulating B/AC-GFs), the atomic in- 
teractions are described by a reactive force field poten- 
tial [36] . ReaxFF is a general bond-order dependent po- 
tential that provides an accurate description of the bond 
breaking and bond formation. Recent simulations on 
a number of hydrocarbon-oxygen systems showed that 
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FIG. 2: (Color online) (a) The spiral B-GO using DFT cal- 
culations for 504 atoms, (b) The spiral B-GO using semiem- 
pirical calculations (PM6) for 1320 atoms. 



ReaxFF rehably provides energies, transition states, re- 
action pathways and reactivity in reasonable agreement 
with ab-initio calculations and experiments. 

For our MD simulations, our rectangular GO nano- 
ribbons contained different numbers of carbon atoms, 
namely Nc = 4800 and 5760. For the GFs we took 
Nc = 5600. For the semi-empirical calculations at the 
PM3 level, we took Nc = 1500. 



III. RESULTS AND DISCUSSION 

A. Infinite one sided hydrogenated graphene 

Figure [1] shows two side views of relaxed infinite (with 
periodic boundary condition applied in both directions 
of the plane) B-GO (a) and AC-GO (b) at r=10K us- 
ing the MD method. In the AC-GO configuration, the 
hydrogens are relatively far from each other and the C- 
H bonds are mutually parallel and perpendicular to the 
z = plane, due to symmetry. In the B-GO the C- 
H bonds are attached to nearest neighbors C-C bond, 
which induces a corrugation in the hydrogen layer with 
C-H bonds that are not fully perpendicular to the z = 
plane, mainly because of the sp^ hybridization but also 
due to the stronger repulsion between te relatively close 
hydrogen atoms. Very similar configurations were ob- 
tained for GF, which are therefore not shown. It is im- 
portant to note that such flat GOs (and one sided GFs) 
are energetically unfavorable [38j. In this study we re- 
port on the structural transformations of these systems 
into stable, non-flat structures. 



B. Ab-initio results: Graphone 

In Fig. [2fa) we show the optimized configuration of B- 
GO using ab-initio calculations with B3LYP/6-31G* and 
504 atoms. Notice that the initial non-relaxed configura- 
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(a) ESP 




(b) HOMO 




(c) LUMO 

FIG. 3: (Color online) Electrostatic potential (a) and the 
HOMO (b) and LUMO (c) for B-GO during rolling. 



tion was a flat sheet. In Fig. El^b) an example of results 
from our semiempirical calculation using PM6 is shown. 

Figure [3l^a) shows the electrostatic surface potential 
(ESP) of B-GO before the completion of the rolling pro- 
cess. The corresponding HOMO and LUMO are shown in 
Figs.[3fb,c) respectively. During rolling both HOMO and 
LUMO are localized at both ends simultaneously which 
is a consequence of the equivalent geometry. The ESP of 
completely rolled B-GO is shown in Fig. |3|^a). The ESP 
indicates that the spiral structure has ionic characteris- 
tics, the red (blue) color indicating positively (negatively) 
charged regions, which stabilizes the spiral geometry even 
more. 

Surprisingly, after complete rolling, the HOMO and 
the LUMO are separated and both localized in only one 
end, as can be seen from Figs. lU^b-c), whereas before 
complete rolling they were localized in both ends simul- 
taneously as already mentioned before (see Figs. [3l^b,c)). 
These effects are related to the differences in electroneg- 
ativity of the simulated atoms (H<C<F). Because of the 
difference in ESP at the two edges, the increment of the 
energy of the LUMO is larger than that of the HOMO, 
which results in a widening of the energy gap from 0.1 eV 
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(b) HOMO 



(c) LUMO 



FIG. 4: (Color online) (a) Electrostatic potential around spi- 
ral graphone and the HOMO (b) and LUMO (c) of the system. 



for the flat sheet to 0.3 eV for the spiral structure. This is 
shown more explicitly in Fig. [51 which reveals an almost 
hnear increase of the gap during the evolution to a spiral. 



C. Ab-initio results: One sided fiuorinated 
graphene 

Very similar results are found for GFs. Figure [6l^a) 
shows the spiral GFs obtained using B3LYP/6-31G*. 
The spiral GF formed using the PM6 method is shown 
in Fig. [6fb). These configurations are minimum energy 
configurations of B-GF. Similarly as for GO, the rolling is 
due to the one-sided, partial (here 50 %) coverage with F 
atoms and corresponding sp^ hybridization. The rolling 
process starts from free boundaries where there is no re- 
sistance against any torque. 



D. Classical molecular dynamics simulation results: 
large size flakes 

Using classical simulations, we first investigated the 
possibility to create hydrogenated CNT from one- 
sided hydrogenated AC-GO nano-ribbons, with periodic 
boundary condition along the zig-zag edge, having a spe- 
cific length equal to L=13.4 A between the two free ends 
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FIG. 5: (Color online) (a) HOMO-LUMO energy diagram 
for four different instances of the rolling process and (b) the 
variation of the energy gap during rolling (for eight steps) for 
the B-GO shown in Fig. EJleft). 
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FIG. 6: (Color online) Spiral B-GF from DFT calculations for 
504 atoms (a). Spiral B-GF from semiempirical calculations 
(PM6) for 1428 atoms (b). 



(in the flat geometry) . After relaxing the system at 10 K 
we found that always a (5,5) CNT is formed independent 
of the dimension of the box in the periodic direction (see 
movie in supplementary material [39]). The used length 
for B-GO, i.e. L=13.4 A, agrees with the perimeter of 
a one-sided hydrogenated (5,5) CNT, i.e. 5^/3ao with 
ao = 1.5 A. The movie in Ref. [3^ shows the time evolu- 
tion for the rolling of an AC-GO which eventually forms 
a (5,5) hydrogenated CNT. This is consistent with pre- 
vious results on the formation of CNT from partially hy- 
drogenated graphene [25]. Note however that the forma- 
tion of these nanotubes is constrained by the choice of 
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(b)HOMO 



(c) LUMO 



FIG. 7: (Color online) ESP(a), HOMO (b) and LUMO (c) 
for the B-GF system shown in Fig. [Gjleft) 



the length of the ribbon, which has to be relatively small 
and should match with the size of an (n,n) tube. 

In the remainder of this section we show the rolling 
effect for much longer nano-ribbons, not leading to tube 
formation, using classical MD simulation. We used gra- 
phone nano-ribbons with free zig-zag edges and periodic 
boundary conditions in the direction parallel to these 
edges (lateral side) and found that the final, minimal en- 
ergy configuration of both GO and GF is a spiral struc- 
ture similar to what we found with DFT. Figures [51^ a, b) 
show two snap shots from the rolled B-GO which are free 
at two longitudinal ends (here set to be zig-zag edges). 
Similar results are found for AC-GOs, see Figs. [Hc,d). 
As expected, the corrugations in B-GO are larger than 
in AC-GO. Thus, unrolling B-GO costs much more en- 
ergy than AC-GO. 

Figure [9] shows the variation of the total energy per 
atom as a function of the rolling process time for AC- 
GO, B-GO, AC-GO* and B-GO*. Clearly B-GO* is the 
most stable having the lowest energy. The energy gap be- 
tween AC-GO (AC-GO*) and B-GO (B-GO*) is about 
0.16 eV/atom. Furthermore it is seen that including the 
vdW energy term lowers slightly the energy by ~ 0.04 
eV/atom due to the attractive van der Waals interac- 
tion. Notice that the ionic characteristic of the rolled 
structure cannot be reproduced by these classical simu- 
lation with AIREBO. Nonetheless, the classically found 
spiral structures are quite stable, which confirms that 
the main driving force for the rolling is the one-sided ab- 



sorption and corresponding sp^ hybridization and not the 
ionic characteristics found in the DFT calculations. 

In Fig. [To] rolled configurations of a system consisting 
of two connected GO parts that are hydrogenated at op- 
posite sides are shown. Hydrogens at x <0 {x >0) are 
bonded to the upper (lower) side of the sheet. Notice 
that the AC-GO and B-GO form different patterns of 
rolled GO with an extra rolled part in the middle. This 
hints to the interesting possibility to engineer the rolled 
shape of GO carefully by selectively hydrogenating parts 
of the initial graphene sheet. 

As an additional, final result of this section we present 
the rolling behavior of a system with one side covered by 
hydrogen and the another side covered by fiuor, which we 
called 'GFH'. Recently, by using ab-initio calculations, a 
negative formation energy was found for this system |40| . 
In Fig. [11] side views are shown for the rolled structure 
of AC-GF (a), AC-GO (b) and AC-GFH (c) after 25 ps 
of molecular dynamics simulations (red, green and blue 
balls are carbon, fiuor and hydrogen, respectively). In all 
three cases we used the ReaxFF potential and the simu- 
lations were done at T=10K, as before. It is clear from 
Fig.[n]that AC-GO (AC-GFH) has the largest (smallest) 
curvature, due to the fact that when one side is covered 
by hydrogen and the other side is covered by fiuor, the 
driven force for rolling is partially balanced, but not com- 
pletely as the fiuor atoms are larger leading to a stronger 
intra-layer repulsion at the fiuor side, so that the system 
still tends to roll, albeit with a much smaller curvature. 



E. Temperature effects 

New structures have to be tested in terms of their sta- 
bility against temperature. In order to study the tem- 
perature effects on the spiral B-GO, after obtaining it at 
10 K (note that in principle spiral GO can also be formed 
at room temperature but we set T=10K in order to min- 
imize thermal fluctuations and to find more rapidly the 
ground state spiral configurations) we raised the temper- 
ature up to 1000 K. Not so surprising in view of the sig- 
nificant gain in binding energy of 0.1 eV/atom ( 1200 K) 
for AC-GO and much more for B-GO (see Fig. 9) in their 
rolled states with respect to the fiat configurations, the 
spiral GOs keep their configuration. Beyond TOOK the 
loops start to deform slightly. However, even at 1000 K 
the system is not un-zipped and the spiral configuration 
is conserved [39| at least within the here used MD simu- 
lation time of 1 ns. For T>1000K the system starts to 
separate into two different loops [39] . From these simula- 
tions we conclude that spiral GO should be rather stable 
at room temperature, facilitating its possible synthesis 
in future experiments. Two movies are provided in the 
supplementary materials to show the rolling of B- 
GO at T=10K and to show the stability of B-GO during 
the heating process up to T=1000 K. For GFs the spiral 
system is somewhat less stable. Close to 800 K the GFs 
start to un-roll and beyond this temperature they start 
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FIG. 8: (Color online) The rolled boat (arm-chair) like graphone nano-ribbons after t=50ps (a,c) and 80 ps (b,d). These 
systems are stable up to at least 1000 K. 
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FIG. 9: (Color online) Variation of the total energy with time 
for both arm-chair and boat-like graphone. 
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FIG. 10: (Color online) The optimized configuration of a 
boat-like (a) and arm-chair like (b) semi-infinite graphone. 
in (a) the A-lattice and in (b) the B-lattice of carbon sites 
were covered by hydrogens. 



to be fractured. Notice that the C-H bonds (C-F bonds) 
in GOs (GFs) are stable even after unrolling, thus no hy- 
drogen release was observed when heating up the rolled 
sheets within our simulation time of 0.1ns. 

We also applied external uniaxial compression along 
the axis of the spiral structures in order to investigate 
the mechanical stability of these structures, which turned 
out to resists against the external pressure much better 
than graphene. Detailed mechanical properties of these 
new spiral hydrocarbons need further investigations. 



IV. CONCLUSIONS 

Free boat like and arm-chair like one- 
sided hydrogenated/fluorinated graphene (gra- 
phone/fluorographene) ribbons/flakes spontaneously 
roll up to form spiral configurations which should be 
quite stable at room temperature. The spiral structures 
exhibit strong mechanical rigidity preventing them 
to unroll. The main driving force behind the spiral 
formation is the energy relaxation associated with 
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FIG. 11: (Color online) Side view of three different spiral 
structures: AC-GF (a), AC-GO(b) and AC-GFH(c) 
after 25 ps of relaxation from the starting flat geom- 
etry. 



the one-sided, asymmetric orientation of sp^ bond- 



ing. Further modestly stabilizing factors are the van 
der Waals attractive interaction between the stacked 
shells and (only in the DFT calculations) the ionic 
interactions between the shells, which appear to be 
polarized across the shell width. Boat-like hydro- 
genated/fluorinated graphene yields more stable, spiral 
systems than arm-chair like graphone. The graphene 
sheet where one side is covered by hydrogen 
and the other side by fluor is unstable and it 
also forms a rolled up structure, albeit with a 
smaller curvature.. The highest occupied and lowest 
unoccupied orbital for spiral GO and GF are localized at 
opposite ends of the system. The energy gap increases 
when the system evolves from the flat to the spiral shape. 
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